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Abstract — We consider the problem of jointly optimizing chan- 
nel pairing, channel-user assignment, and power allocation, to 
maximize the weighted sum-rate, in a single-relay cooperative 
system with multiple channels and multiple users. Common relay- 
ing strategies are considered, and transmission power constraints 
are imposed on both individual transmitters and the aggregate 
over all transmitters. The joint optimization problem naturally 
leads to a mixed-integer program. Despite the general expectation 
that such problems are intractable, we construct an efficient 
algorithm to find an optimal solution, which incurs computational 
complexity that is polynomial in the number of channels and the 
number of users. We further demonstrate through numerical 
experiments that the jointly optimal solution can significantly 
improve system performance over its suboptimal alternatives. 



I. Introduction 

We consider the problem of resource assignment for multi- 
channel multi-user communication through relaying. The prob- 
lem typically arises in cellular communication or wireless local 
area networks, through either dedicated relay stations or users 
temporarily serving as relay nodes. In traditional narrow-band 
cooperative relaying systems, the relay retransmits a processed 
version of the received signal over the same frequency channel. 
In contrast, when multiple frequency channels are available, 
the relay can exploit the additional frequency dimension, to 
process incoming signals adaptively based on the diversity in 
channel strength. 

In narrow-band cooperative relaying systems, the relay 
retransmits a processed version of the received signal over the 
same frequency channel. In contrast, when multiple frequency 
channels are available, the relay can exploit the additional 
frequency dimension, to process incoming signals adaptively 
based on the diversity in channel strength. Channel pair- 
ing, which devises a matching of incoming and outgoing 
subcarriers in OFDM-based relaying, was proposed indepen- 
dently in Q and for single-user relayingfj In a multi- 
user communication environment, both incoming and outgoing 
channels at the relay are shared among all users. A cru- 
cial problem is to determine the assignment of a subset of 
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1 Since a vast majority of multi-channel relaying systems in the literature 
are based on OFDM, we use it as an illustrative example in this work, so that 
the terms "channel" and "subcarrier" are synonymous. 



incoming-outgoing channel pairs to each user, which we term 
channel-user assignment. Since the channel condition can vary 
drastically for different users, and over the same incoming 
and outgoing channels, judicious channel-user assignment and 
channel pairing can potentially lead to significant improvement 
in spectral efficiency. Together with power allocation over 
multiple channels at the transmitters, essential for performance 
optimization, these are three main resource assignment prob- 
lems in multi-channel multi-user relaying. 

There is strong correlation among channel pairing, channel - 
user assignment, and power allocation. Joint consideration of 
these three problems is required to achieve optimal system 
performance. However, the combinatorial nature of channel 
pairing and assignment generally leads to a mixed-integer 
programming problem, whose solution often bears prohibitive 
computational complexity and renders the problem intractable. 
As a result, previous attempts to optimize the performance 
of multi-channel multi-user relaying systems through resource 
allocation often consider only a subset of these three problems 
- lfl6l , or adopt suboptimal approaches |fl7l - ||23l . 

In this work, we consider all three resource assignment 
problems in a dual-hop multi-channel relaying network for 
multi-user communication through a single relay, under several 
common relaying strategies. We show that there is an effi- 
cient method to jointly optimize channel pairing, channel-user 
assignment, and power allocation in such general dual-hop 
relaying networks. The proposed solution framework is built 
upon continuous relaxation and Lagrange dual minimization. 
Although this approach is often applied to integer program- 
ming problems l24l . it generally provides only heuristic or 
approximate solutions. However, by exploring the rich struc- 
ture in our problem, we show that judicious reformulation 
and choices of the optimization trajectory can preserve both 
the binary constraints and the strong duality property of the 
continuous version, thus enabling a jointly optimal solution. 

Through reformulation, we transform the core of the orig- 
inal problem into a special incidence of the class of three- 
dimensional assignment problems, which is NP hard in gen- 
eral but has polynomial-time solutions - in terms of the 
number of channels and users - for our specific setting of 
channel pairing and channel-user assignment. For the often 
studied conventional decode-and-forward (DF) relaying with 
a maximum weighted sum-rate objective, we further propose 
a divide-and-conquer algorithm for dual minimization, which 
guarantees that convergence to an optimal solution requires 
only a polynomial number of iterations in the number of 
channels. This ensures the scalability of the proposed solution 
to large multi-channel systems. 
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Our proposed solution is applicable to a wide range of sys- 
tem configurations, accommodating both total and individual 
power constraints, and allowing direct source-destination links 
in relaying. We show that it can be modified to work with var- 
ious relaying strategies in addition to DF, including variants of 
compress-and-forward (CF) and amplify-and-forward (AF). It 
also accommodates general concave utility functions. Through 
simulation and numerical comparison, we further illustrate that 
there is often a large performance gap between the jointly 
optimal solution and the suboptimal alternatives. 

The rest of the paper is organized as follows. We first 
provide a literature review of the related work in Section [EI] 
In Section [Till we discuss the system model and formulate 
the joint optimization problem. For weighted sum-rate max- 
imization with DF relaying, we describe our framework of 
finding the optimal solution with polynomial complexity in 
Section[lV] Extension to other relaying strategies are explained 
in Section [V] Numerical studies are presented in Section |VI] 
and conclusions are given in Section IVIII 

II. Related Work 

Most existing works on optimizing resource allocation for 
multi-channel relaying systems consider a subset of the three 
aforementioned problems. After channel pairing was proposed 
in 12 and Q, its optimization has been considered in several 
studies. In the absence of the direct source-destination link, J4] 
showed that the sorted-SNR channel pairing scheme, which 
matches the incoming and outgoing subcarriers according 
to the sorted order of their SNRs for some given power 
allocation, is sum-rate optimal for a single-user AF relaying 
OFDM system. When the direct source-destination link is 
available, a low complexity optimal channel pairing scheme 
was established in J5) for AF relaying. In addition, it was 
shown that channel pairing is optimal among all unitary linear 
processing at the relay under a fixed gain power assumptions. 
However, none of these works considered optimizing power 
allocation. Channel-user assignment in multi-user relaying 
networks, under given power allocation, was considered in 
0, where the authors sought an optimal channel-user and 
channel-relay assignment to maximize the uplink data rate for 
AF and DF relaying with multiple relays. For a multi-channel 
network with multiple sources, single AF relay, and single 
destination, Q studied the problem of channel pairing and 
channel-user assignment. It maximizes the sum received SNR, 
assuming that the power allocation is given. A suboptimal 
solution is proposed for distributed implementation using game 
theory. Finally, the problem of optimal power allocation for 
OFDM relaying in specific relay network setups was studied 
in numerous works for different relay strategies and power 
constraints, see for example J8], J9], ifTO) . 

Jointly optimizing channel pairing and power allocation for 
single-user relaying was considered in several studies. Without 
the direct source-destination link, ifTTI and |[T2l considered this 
problem for dual-hop DF relaying in an OFDM system for 
total power and individual power constraints, respectively. It 
was shown that joint subcarrier pairing and power allocation 
are separable for sum-rate optimization. This separation was 



also established in the general multi-hop case in |[T3l . for 
both AF and DF relaying, and under either total power or 
individual power constraints. With consideration for the direct 
source-destination link, the authors of |[l4l and lfT31 studied 
joint subcarrier pairing and power allocation in a single-user 
OFDM system, for AF and DF relaying respectively. The 
joint optimization problems were formulated as mixed-integer 
programs and solved in the Lagrange dual domain. Although 
strict optimality was not established, the proposed solutions 
were shown to be asymptotically optimal as the number of 
subcarriers approaches infinity, based on the frequency-domain 
virtual time-sharing argument ll25l . For relay-assisted multi- 
user scenarios, joint optimization of channel-user assignment 
and power allocation was considered in |[T6) for communica- 
tion between a base station and users who have the ability to 
relay information for each other. Based on the same virtual 
time-sharing argument ll25l . asymptotically optimal solution 
was provided for network utility maximization. 

The problem is especially challenging when channel pair- 
ing, channel-user assignment, and power allocation need to 
be optimized jointly in relay-assisted multi-user scenarios. 
Existing work to tackle it has been scarce. In IfTTI . such joint 
optimization was considered for cooperation among users in 
uplink communication, accounting for the splitting of band- 
width at a user that needs to simultaneously transmit its own 
data and relay for others. The proposed problem was NP hard 
and a suboptimal heuristic algorithm was constructed. The 
authors of IT8l studied this problem for a single relay using 
DF without the direct source-destination link. Under a total 
power constraint, they showed that, for sum-rate maximization, 
it is optimal to separately design channel-user assignment, 
channel pairing, and power allocation. However, this approach 
is suboptimal for the general case when the direct link is 
available, when the user weights are non-uniform, or when 
individual power constraints are considered. In comparison, 
we consider more general relaying strategies that use the direct 
source-destination link, so that the simple pairing scheme 
based on sorted channel gain is no longer optimal. Further- 
more, our proposed approach accommodates individual power 
constraints in addition to total power constraints, relaying 
strategies other than DF, and other optimization objectives. 
In Section IVII we further illustrate with numerical data that 
there is a large performance gap between such a separate 
optimization approach and the jointly optimal solution. 

There are also other studies on resource allocation in multi- 
channel relaying systems, with different system models from 
the one presented in this paper (for example, |fl9l , |f20l . lETTl . 
JUL l23lD - Due to the significant complexity in these system 
models, no general optimal solution has been found. Rather, 
suboptimal algorithms are proposed with an aim to support 
satisfactory system performance. In contrast, in this work 
we tackle the problem of joint resource optimization in a 
simpler, single-relay system with multiple users, proposing a 
provably optimal solution with a formal proof for polynomial- 
time complexity. Some preliminary results of this study have 
appeared in JT]. This version contains substantial extensions, 
adding detailed solutions on how to accommodate alternate 
power constraints and performance objectives, and presenting 
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Fig. 1. Illustration of dual-hop multi-channel relaying. 

new derivations, proofs, and numerical results. 

III. System Model and Problem Formulation 

We consider the scenario where a source communicates with 
K users via a single relay as illustrated in Fig. Q] The available 
radio spectrum is divided into N equal bandwidth channels, 
accessible by all nodes. We focus on the downlink in our 
analysis in this paper, but the proposed solution framework 
can be adopted for the uplink by swapping the roles of the 
source and the users. 

We denote by hf, h rk , and hf k the state of channel i, for 
1 < i < N, over the first hop between the source and the relay, 
over the second hop between the relay and user k, and over the 
direct link between the source and user k, respectively. The 
additive noise on a channel at the relay and user k are modeled 
as i.i.d. zero-mean Gaussian random variables with variances 
of and <t2, respectively. The channel state is assumed to be 
available at both the source and the relay, which enables them 
to dynamically assign channels and allocate power according 
to channel conditions. 

A. Channel Assignment 

The relay transmits a processed version of the incoming data 
to its intended user using a specific relay strategy. The relay 
also conducts channel pairing and channel-user assignment. 
Channel pairing refers to a one-to-one mapping between the 
incoming channels and outgoing channels at the relay. Through 
channel-user assignment, on the other hand, a subset of 
incoming-outgoing channels is assigned to each user. Clearly, 
channel-pairing choices are closely connected with how the 
channels are assigned to the users. We term the joint decision 
on channel pairing and channel-user assignment the channel 
assignment problem. As the different channels exhibit various 
quality, judicious channel assignment can potentially lead to 
significant improvement in spectral efficiency. 

We say a path V(m, n, k) is selected, if first-hop channel m 
is paired with second-hop channel n, and the pair of channels 
(m, n) is assigned to user k. We define indicator functions 
4>mnk for channel assignment as follows: 



1, if V(m,n, k) is selected, 
0, otherwise . 



(1) 



There is a one-to-one mapping between first-hop channels 
and second-hop channels. Furthermore, we require that each 
channel pair be assigned to only one user, but a user may be 



assigned multiple channel pairs. Hence <j) mn k is constrained 
by 



JV K 



N K 



y <t>mnk = i, Vm, ^2 E ^ mnk = ^ n - ( 2 ) 



n=l k=l 



B. Power Allocation 



m—1 k—1 



Along any path V(m, n, k), the source and relay transmis- 
sion powers are denoted by P r s nnk and P£ mk , respectively. We 
consider both individual power constraints, 



N N 



K 



N N 



K 



EEE^ 

n—1 m—1 k—1 



< P, 



EEE p ™ 



■mnk 



< p r 



=1 m=l fc=l 



and the total power constraint, 

N N K 
n—1 m—1 k—1 



pr \ <■ P 
£ mnk) — r t\ 



(3) 



(4) 



where P s , P r , and P t are the maximum allowed transmission 
power by the source, the relay, and the combined source 
and relay, respectively. This is a general representation of 
the power limitations imposed on the system including, e.g., 
hardware constraint, legal or regulatory requirement, or energy 
conservation. Note that this general representation can be 
easily tailored to also specify systems with only individual 
power constraints, or only total power constraint, by setting 
one or more of P s , P r , and P t to sufficiently large values. 

Each constraint above is either inactive (i.e., at optimality it 
is satisfied with strict inequality) or active (i.e., at optimality 
it is satisfied with equality). We consider the case where all 
active constraints are strictly active, i.e., if the problem is 
modified by changing the power limits by small amounts, at 
optimality the constraints remain active. This is without loss 
of generality, since any constraint that is active but not strictly 
active can be made inactive, by increasing the power limit by 
a small amount, without altering the problem solution. 
Define p mnk 

(Ps,Pr,Pt)- 



/ ps pr ps 

V mnk ' mnk ' mnk 



Pmnk) ^d 7T = 



C. Relaying Strategy 

We initially focus on DF relaying but will later show how 
the proposed method can be applied to other relaying schemes, 
such as AF and CF. We consider a general case where, apart 
from the relay path, the direct links are available between the 
source and users. In this case, the signals received from the 
relay path and the direct link can be combined to improve the 
decoding performance. In DF, each transmission time frame 
is divided into two equal slots. In the first slot, the source 
transmits an information block on each channel, which is 
received by both the relay and the intended user. In the second 
slot, the relay attempts to decode the received message from 
each incoming channel (first hop), and forwards a version of 
the decoded message on an outgoing channel (second hop) 
to the intended user. The intended user collects the received 
signals in both time slots, applies maximum ratio combining, 
and decodes the message. 
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Consider the conventional repetition-coding based DF re- 
laying ||26l , |27l , where the relay is required to fully decode 
the incoming message, re-encode it with repetition coding, 
and forward it to the intended user. The maximum achievable 
source-destination rate on path V(m,n, k) is given by ll27l 

R(m, n, k) = - min{log(l + a m P^ nk ), 



log(l 



CmkPmnk 



(5) 



where a„ 



~Zi — , °nk 



\K 



\h sh \ 2 

and c rnk = are 



normalized channel power gains against the noise variance at 
the relay and user k, and the base of logarithm is 2. 

D. Optimization Objective 

Various rate-utility functions can be used as objectives. For 
convenience of illustration, in this paper we focus on the 
weighted sum-rate. Denoting by w k the relative weight for 
user k, such that Y^=i Wk ~ 1> we formulae the problem of 
weighted sum-rate maximization as 

K N N 

max w k y^ (f> mnk R{m,n,k) (6) 

k—1 ri— 1 rn— 1 

s .t. da,®,©, 

(f>mnk G {0, 1}, \fm,n,k (7) 
P^ nk >0, P mnk >0, Vro,n,fc, (8) 

where * = [(/) m nk}NxNxK, P s = [P^, nk ]NxNxK, and P r = 
[PmnklNxNxK- Given the relative weights and the channel 
gains on each path V(m,n, k), the optimization problem © 
finds the jointly optimal solution of channel pairing, channel- 
user assignment, and power allocation by optimizing <&, P s , 
and P r . 

IV. Weighted Sum-Rate Maximization for 

MULTI- CHANNEL DF 

The optimization in ([6) is a mixed-integer programming 
problem, which in general has intractable complexity due to 
its combinatorial nature. However, in this section, we present 
a method to find an optimal solution with computational 
complexity growing only polynomially with the number of 
channels and users. 

A. Convex Reformulation via Continuous Relaxation 

The proposed approach is built on the reformulation of (|6} 
into a convex optimization problem with a real-valued <£> and 
strong Lagrange duality. We later show that the reformulated 
problem is optimized by a binary 3? = <fr. 

We first substitute 

pr 

pr mnk /q\ 

r mnk — ~ 



ps 

ps mnk 



4^mnk 



and 



Vmnk 



into the objective of ©. This does not change the original opti- 
mization problem, since if <j> mn k = 1, then (0 is trivially true; 
and if <fi mnk = 0, then by l'Hopital's rule, cf> mnk R(rn,n,k) 
remains zero before and after the substitution. Indeed, it obvi- 
ously preserves the optimality of power allocation to enforce 



We then relax the binary constraint on 3? by defining a 
continuous version of 4> mnk , denoted by (j> mn k, which may 
take any value in the interval [0,1]. Then, the reformulated 
version of the optimization problem © can be written as 

; mm{log(l + a m 5^), 



-—( 

, 2 



J mnk 



P, 



log(l + c mfcT 



mnk 



pr 



s.t. ^ 4>mnk = l,Vm, ^ q 

n,k m.k 

< 4> m nk < 1, Vm,n,k, 

©,©,©. 



+ b nk ~f^)} (10) 

mnk tymnk 

l,Vn, 



(11) 



(12) 



The objective function ( fTOb is concave in (<&,P s ,P r ), 
since cf> mnk log(l + a m j^) and 4> mnk log(l + c m kj^; + 

pr 

b nk ■r m7lk ) are the perspectives of the concave functions 



log(l_ 



" a ™P mnk ) and log(l + c mk P° 



b nkP mnk ), respec- 



' \ |— | ■ ■ - 1 1 L 1 b r\i / ' \ ■ ■ ■ I l i t ' ri ■ I tli I on / a 

tiveljQ It is also noted that the minimum of two concave 
functions is a concave function. Furthermore, since all the 
constraints are affine, and there are obvious feasible points, 
Slater's condition is satisfied l28l . Hence, the convex opti- 
mization problem ([Tol l has zero duality gap, suggesting that a 
globally optimal solution can be found in the Lagrange dual 
domain. 

Using continuous relaxation on integer programming prob- 
lems is not a new technique (24|. However, doing so typically 
leads only to heuristics or approximations. Clearly, solving a 
maximization problem with relaxed constraints generally gives 
only an upper bound to the original problem. In particular, all 
global optima for ([Tol l do not necessarily give a binary <i>, 
which is required for (|6). However, we next show that, in 
the problem under consideration, indeed there always exists a 
globally optimal solution to ( fTOb consisting of a binary <&, and 
the proposed approach ensures that such an optimal solution 
is found in polynomial time. 

B. Power Allocation via Maximization of Lagrange Function 
over P s and P r 

Consider the Lagrange function for (TTOl . 
£($,P s ,P r ,A) = ^4>mnk min{log(l + a n 1 



Vrnnk 



log(l + c„ 



+b nk 5^)}-(\ 3 +A t ) J2 k 



k^r 



{K + At) J2 Plunk + ATT 2 



(13) 



ps _ pr 

mnk mnk 



for all (m, n, k) such that 



0. 



where A = (A s , A r , At) is the vector of Lagrange multipliers 
associated with the power constraints (O and (0J. The dual 
function is therefore 

g{\) = max £(*, P s , P' , A) (14) 

ps pr 

2 The perspective of function / : R™ — > R is defined as g(x,t) = tf(x/t), 
with domain {(x,t)\x/t G dom/, t > 0}. The perspective operation 
preserves concavity 1281 . Here we include 4> mn k = in the domain of 
the perspectives. It is easy to see that they remain concave. 
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s.t. CED, CE3, ©. 

The above maximization of the Lagrange function can be 
carried out by first optimizing the power allocation given fixed 
<I>. The KKT conditions suggest that the maximization of 
( TBi i over P s and P r can be decomposed into N x N x K 
independent subproblems to find the optimal P^ nk and P m * nk - 



max 

p', >o,pr 



>o 



mnk \ 



t>mnk } P m nk ' Pmnk I A) 



(15) 



where £ mnfe (> m „ fe , P^ nfc , P^ nfc , A) is the part of 
£(<&, P s , P r , A) that concerns only the path Vim, n, k). 

It can be shown that the solution to (fT~5b has the following 
form. The derivation details are given in Appendix [A] Note 
that, since P m * nk and P m * nk depends on 4> mn k in an obvious 
way, we simply present them as functions of A for the rest of 
this section. 



1 



a(X s +X t ) 

, if a m 



0mnk, 



Pi 

arg max £ 
L(p=,p-)e{ Pl ,p 2 } 



^ C-mk 

mnk 



and 



rank ; 
c m fc 



< 



if a n 



X s +Xt ' x r +x t 
k ,P s ,P r ,X), o.w. 



(16) 



-'inn 



where a = 2 In 2, [x) + = max{x, 0}, pi 

W k Dn 



{b nk (X 3 +X t ) + (a 

Wk 



a(X s +X t ) 



1 

Cmk 



c mk )(X r +X t )) 
firnnk, 



'Jmnki 



and P2 = 



C. Channel Assignment via Maximization of Lagrange Func- 
tion over $ 

To maximize the Lagrange function over <&, we define 
1 



^-mnfc(A) 



"mnk 



-£ 



mnk \(ymnk j 



Pmnk W J Pmnk (A) ! A) 



(17) 



Note that A mn fc(A) is independent of 0„mfc because of the 
multiplication form of ( fTST l by (pmnk- Then, (TBI can be 
determined by the following optimization problem over 4>: 



max <p 

in ilk - 1 /// // A: (A) 



(18) 



.t. CD,©- 



To proceed, we present the following lemma on the decom- 
position of 3>. 

Lemma 1: Any matrix 3? = [4>mnk\NxNxK with < 
4>mnk < 1 and satisfying (fTTT i can be decomposed into one 
matrix X = [:r mn ]jv x jv and MN vectors y m " = [y k nn ]i X K, 
such that mnfe = cc mn yj!"\VTO,n, fc, with < x m „ < 1 
and < < 1, satisfying J2 n x mrl = l,Vm, £) m ^mn = 
1, Vn, and ^ fc y™" = l,Vm, n. Furthermore, any such matrix 
X and vectors y mn uniquely determines a matrix 3? that is 
given by <f> mnk = x mn y k nn and satisfies ([[]}■ 

Proof: The proof is provided in Appendix [B] ■ 

Note that, even though the above decomposition can also be 
applied to a binary $ as a trivial special case of Lemma Q] we 



require the general form of this lemma to deal with continuous 
values in <&, X, and y mn . In particular, the mapping from <J> 
to (X, {y mn }) is one-to-many, which is quite different from 
the binary case. 

Lemma [T] implies that any optimization over (X, y" m ) 
also optimizes 3? for the same objective. This allows us to 
replace, in problem (11 8t . <f> mnk with x mn ?/™™. Furthermore, 
the constant terms can be dropped from (ITFt . Hence, we can 
equivalently seek solutions to the following problem 



max V" Xmn V" Z/r™ Annfc (A) 
lJ ' m,n k 



(19) 



S.t. ^ x mn = 1, Vm, 2J ^mn = 1, Vn, < X mn < 1, Vm, 1 

n m 

(20) 

J2vr = ^m,n, 0<yr<l,Vrn,n,k. (21) 
fe 

The following two-stage solution is sufficient. First, the inner- 
sum term is maximized over y k nn for each (m, n) pair, i.e., 



(22) 



y i 

fc 

s.t. d2lT ). 

An optimal solution to d22b is readily obtained as 

y mn* = f 1 > if fc = ar S max l<Z<K A TO „;(A) (23) 

[ 0, otherwise 

In the above maximization, arbitrary tie-breaking can be 
performed if necessary. Next, inserting ^^(A) into ( fl9l ). we 
have the linear optimization problem 



max x mn A mn ( A) 

m,n 



(24) 



It is well known that there always exists an optimal solution to 
( l24l ) that is binary |24] Chapter 3]. An intuitive explanation is 
the following. Since d24b is a linear program with a bounded 
objective, an optimal solution can be found at the vertices 
of the feasible region. Furthermore, since X is a doubly 
stochastic matrix, it is a convex combination of permutation 
matrices. One of these vertex permutation matrices is an 
optimal solution to ([24}, so at least one optimal X is binary. 
Then, to find a binary optimal X, d24b is a two-dimensional 
assignment problem. Efficient algorithms, such as the Hungar- 
ian Algorithm 1291 . exist to produce an optimal solution with 
computational complexity being polynomial in N. 
Finally, the optimal <f> mnk given A is 



Cn*(A) = aSm(A)ylT*(A) 



(25) 



Since binary x mn (X) and j/™ n *(A) are computed following 
the above procedure, 4>* mnk {\) is also binary. This shows 
that there exists at least one binary optimal solution to the 
maximization in ( TT~8b . 

Intuitively, the globally optimal solution described above 
suggests a pairing between the input and output channels at 
the relay, and if channels m and n are paired, they are assigned 
to a single user k, whose associated A mnk (X) is the greatest 
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among all users. Note that such an interpretation might lead us 
to conclude that we could have forgone continuous relaxation 
from the very beginning and focused only on a binary <1>. 
However, we would still have required the continuous <& to 
construct a convex optimization problem, whose strong duality 
property provides the optimality of the proposed approach. The 
optimality of (X, {y mn }) taking binary values is implied only 
through the above derivation. 

Interestingly, the original optimization problem (ITS) with 
a binary matrix $ is a special case of the axillary three- 
dimensional assignment problems 11301 . It is well known that 
the general form of this family of problems is NP hard and 
cannot be solved by continuous relaxation on unlike the 
two-dimensional assignment problem in (l24l . In our case, the 
special structure of 3> expressed in ill U . namely the absence 
of a constraint on per-user resource allocation, makes possible 
the availability of an efficient solution to (1181 1. 

It is also worth noting that, given any A, there may exist 
non-integer optimal solutions to (fT~8T >. For example, when the 
maximal value of A mn k in ( 1221 is achieved by multiple users 
having the same channel gains, there is an infinite number 
of optimal y mn , leading to non-integer optimal solutions 
for (fimnk- However, the procedure above finds only one of 
the optimal solutions in binary form, which is sufficient for 
computing the dual function. 

D. Dual Minimization: Baseline Subgradient Approach 

The previous subsection provides a way to find the Lagrange 
dual g(X) for any Lagrange multiplier vector A. Next, the 
standard approach calls for minimizing the dual function: 



min g(X) 
s.t. A X 0. 



(26) 



This can be solved using the subgradient method PP . It is 
easy to verify that a subgradient at the point A is given by 



0(A) =7T- Y, P*mnkW . 



(27) 



where P^* nk {X) and P^nkW are computed based on ([Tfjl i 



and tfnnkW found usin 8 

For completeness, we first summarize the standard sub- 
gradient updating algorithm for solving the dual problem in 
the following. We will present a modified dual minimization 
algorithm in Section IIV-FI which is guaranteed to converge in 
polynomial time. 

1) Initialize A (0) . 

2) Given A (i) , obtain the optimal values of P^ nk (X^), 
P^A(<>), and 0* (AW). 

3) Update A through A (i+1) = [A m - 0(A w )z/°] + where 
i/W is the step size at the kh iteration. 

4) Let 1 = 1+1; repeat from Step 2) until the convergence 
of min; g(X"'), 

Several step-size rules have been proven to guarantee con- 
vergence under some general conditions ll3T1ll32l . For exam- 
ple, using a constant step size v, i.e., = v, or using a 
constant step length v, i.e., = ^/||0(A ( -'' > )||2, leads to an 



objective within a given neighborhood of a global optimum; 
while using the non-summable, square-summable rule leads 
to asymptotic convergence to a global optimum. Furthermore, 
one may satisfy any constraints on X^ within a convex region 
by projecting A^ onto the region. This is the general projected 
subgradient method, which does not reduce the speed of 
convergence l32l . For example, Step 3 above ensures that 
A''' y 0, and we will further consider projection onto convex 
regions TZ\ and IZ2 in Section IIV-FI 

E. Primal Optimality 

With standard subgradient updating, the dual optimal A* is 
obtained, from which we compute the channel assignment and 
power allocation matrices (<£*, P s *, P"), where $* = <fr*. 
Since the optimization problem ( fTOb is a convex program that 
satisfies Slater's condition, it has zero duality gap. Denote 
by f*(ir) the maximal value of the objective in (fTOb - Then 
f*(n) = g(X*), and furthermore it is concave in n. We 
consider systems that have the following strictly diminishing 
rate-power relation: 

Assumption 1: /*(tt) is strictly concave in any strictly 
active power constraint P x <G {P s , P r , P*}. 

In other words, as the data rate increases, each unit of 
increment requires more and more marginal power. With a 
strictly concave R(m, n, k) in terms of p mnk , this assumption 
holds when either there is no tie-breaking in (l23~t or (l24l or 
there is tie-breaking that is due to users or paths having the 
same weights or channel gains^. 

Proposition 1: Under Assumption Q] (<£*, P s * ; P") is a 
globally optimal solution to the original problem ©. 

Proof: Since P s * and P s * are uniquely determined by 
A* and f we need only to focus on 3?*. 

For any inactive constraint P x , we have X* = and the 
subgradients of g(X) in the direction of X x are all positive. 
Hence any 3>* is feasible with respect to P x . 

For any strictly active constraint P x , we have A* > 0. 
Furthermore, 



/*(tt) =£(**, P S *,P", A*) 



(28) 



m,n,k 



m.n,k 



Given Assumption Q] the above is possible only when all 
strictly active constraints are satisfied with equality. Therefore, 
any is feasible with respect to P x and the complementary 
slackness condition is satisfied. Hence, ($*,P s *,P r *) is a 
globally optimal solution to (fTOb . 

Furthermore, since ([T(7t is a constraint-relaxed version of 
©, <&* gives an upper bound to the objective of ([6j. Finally, 
since 3>* satisfies the binary constraints in © at each iteration 
of the subgradient algorithm, it satisfies all constraints in ©. 
Therefore, it is a globally optimal solution to ©. ■ 

We point out that using conventional convex optimization 
software packages directly on the relaxed problem (TTOt is not 

3 However, we cannot rule out the possibility of a case where other forms of 
tie-breaking in j23j or J241 might create linear segments in /*(7r), although 
in all simulation tests with arbitrary parameters, we have not produced a case 
where this assumption fails. 
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sufficient to solve ©. This is because there is no guarantee 
that they will return a binary <&*, and furthermore due to 
complicated three-dimensional dependencies among <j> mn k, 
there is no readily available method to transform a fractional 
<J>* to the desired binary solution. 

F. Dual Minimization: Divide-and-Conquer Algorithm with 
Polynomial Complexity 

The standard subgradient method produces a global op- 
timum, but its computational complexity is not generally 
known. Previous studies have provided asymptotic bounds or 
conjectures on its efficiency through computational experience. 
In general, the number of iterations in subgradient updating 
depends on the step-size rule, the distance between the initial 
solution and the optimal solution, and the 2-norm of the 
subgradients 1311, l3"2l. 

Next, we propose a new dual minimization algorithm that 
guarantees convergence with polynomial complexity in N and 
K, to a global optimum for our optimization problem. It 
uses a divide-and-conquer approach, by grouping the possible 
locations of A* into two regions and applying projected 
subgradient updating constrained within either. It ensures that 
in each region, our choice of the initial A*- ' and subsequent 
subgradient updating lead to convergence in polynomial time. 

We first define the following two overlapping convex re- 
gions in terms of A: 



A 



A : A s + A t > 



min Wk min{ min a m , min c m k} 

{k:w k >0} {m:a m >0} {m,k:c mk >0} 

4a(maxa m min{P s , P t ] + 1) 
min a m 

n 2 ^\x:x s + x t + {m - a - >0 h } (X r + X t )> 

max bnk 

n,k 

min Wk 

{k:w k >0} 



a y o 



, a y o 



a(min{P s , P t } 



mm a Tl 

{m:a m >0} 



These are two possible regions where A* resides, which 
depends on whether there exists at least one chosen path with 
non-zero direct-link channel gain c m k- This is formalized in 
the following lemma. Its proof is given in l33l . 

Lemma 2: If there exists some (m, n, k) such that 4>* nnk = 
1 and Cmk > 0, then A* £ TZ\. Otherwise, A* £ TZ 2 . 

The proposed divide-and-conquer dual minimization 
(DCDM) algorithm considers both possible regions for A*. 
It first creates the two conditions in Lemma [2] by artificially 
setting direct-link channel gains to zero. It then applies the 
projected subgradient algorithm on Tvli and TZi separately, and 
chooses the better solution between these two. The algorithm 
is formally detailed in Algorithm Q] and its optimality and 
complexity are given in Propositions [2] and [3] respectively. 
Note that one cannot use Lemma [2] to determine, before the 
optimal channel assignment matrix is chosen, which region 
A* is in. This necessitates the comparison step in the DCDM 
algorithm. 



Algorithm 1 Divide-and-Conquer Dual Minimization 
(DCDM) 

if there exists some m and k such that c m k > then 
AJ = output of subgradient updating algorithm with 
projection onto A^ £ IZi 

Set c mk = for all 1 < m < N and 1 < k < N 

Aj = output of subgradient updating algorithm with 

projection onto A^ £ IZ2 

return argmin <?(A) 

Ae{A*,A*}' 

else 

A* = output of subgradient updating algorithm with 
projection onto A*-'' £ IZ2 
return A* 
end if 



Proposition 2: With DCDM, the computed channel assign- 
ment and power allocation matrices (<&*, P 3 * , P"), where 
$* _ <j>*^ j s a globally optimal solution to the original 
problem ©. 

Proof: Suppose there exists some (m, n, k) such that 
't'mnk = 1 an d c «ifc - > 0. Then Lemma [2] shows that 
A* £ Hi, Therefore, by Proposition [T] AJ obtained by 
subgradient updating projected onto TZi is an optimal solution. 
Furthermore, setting c mk = for all 1 < m < N and 
1 < k < N only reduces R(m, n, k) for all paths, so that 
subsequently minimizing the Lagrange dual yields an inferior 
solution. Therefore, argmin A6 | A » ^»\fl l (A) = A* is returned 
by DCDM. 



Suppose c 



/ ) I k 



for all (m, n, k) such that 



k = 1. 

' i.e., all chosen paths have zero direct-link channel gain. Then, 
setting c m k = for all 1 < m < N and 1 < k < N only 
reduces R(m, n, fc) for the non-chosen paths. Subsequently 
minimizing the Lagrange dual yields the same solution as 
before changing c m k- Furthermore, Lemma |2] shows that this 
optimal solution is in TZ 2 - Hence, \^ obtained by subgradient 
updating projected onto 1Z 2 is an optimal solution. In this case, 
argmin Ag | A » A ,j <?(A) = A2 is returned by DCDM. ■ 

The polynomial computational complexity of DCDM is 
stated in Proposition [3] Its proof requires the following lem- 
mas, which give upper bounds on j|A*|| 2 and ||0(A (i) )|| 2 , 
where || ■ H2 denotes the 2-norm. 

Lemma 3: At global optimum, || A*|| 2 is upper bounded by 
A ma , = 0{N 2 ). 

Proof: The proof is provided in Appendix [D] ■ 

Lemma 4: At every step of subgradient updating in the 
DCDM algorithm, ||0(A W )||2 is upper bounded by 6 max = 
0(N 2 ). 

Proof: The proof is provided in Appendix [E] ■ 
Proposition 3: To achieve a weighted sum-rate within an 
arbitrary e > neighborhood of the optimum <?(A*), using ei- 
ther a constant step size or a constant step length in subgradient 
updating, the DCDM algorithm has polynomial computational 
complexity in N and K. 

Proof: At each iteration of the standard subgradient 
updating algorithm, the procedures described in Sections IIV-BI 
and IIV-CI are employed. This has computational complexity 
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polynomial in TV and K. Therefore, it remains to show that 
the total number of iterations is not more than polynomial in 
N or K. 

For either case of projecting onto 1Z\ or TZ.2, one may choose 
an initial A*- -* such that the distance between A*- -* and A* is 
upper bounded by X ma x- Then, it can be shown that, at the Zth 
iteration, the distance between the current best objective to the 
optimum objective g(X ) is upper bounded, by 2 vi m ° x 
if a constant step size is used (i.e., i/W = v), or by 



A; 



2iA 

COM 



if a constant step length is used (i.e., z/W 



is I \\0{\ (l) ) || 2 ) E]||32). For the former and latter bounds, if we 
set v = e/9 max and v = e/9 rnax respectively, both are upper 



J max 

bounded by e when I > \ max 0^ nax /e' i = 0(iV 4 ). Hence, the 
number of required iterations until convergence, for either of 
the two projected subgradient updating procedures in DCDM, 
is polynomial in N and independent of K. ■ 
Note that using the non-summable, square-summable step- 
size rule in the early iterations of subgradient updating, often 
leads to faster movement toward a global optimum than using 
a constant step size or a constant step length. This is due to its 
larger step sizes when I is small. However, such a step-size rule 
does not guarantee polynomial convergence time0. Therefore, 
one may start with the non-summable, square-summable rule, 
and then switch to one of the constant-step rules when the 
step size or step length is sufficiently near the prescribed 
value in Proposition [3] This would reduce the convergence 
time in practice while preserving the guarantee of polynomial 
complexity. 

V. Extensions to General Relaying Strategies 

For any relaying strategy in which data sent through 
different communication paths Vim, n, k) are independent 
and the achievable rates R(m, n, k) is a concave function 
in transmission powers (Pm„k,Pmnk)' tne P ro P ose d solution 
approach gives jointly optimal channel assignment and power 
allocation for weighted sum-rate maximization. To see this, we 
first note that any concave rate function would lead to convex 
programming for the relaxed and reformulated problem, which 
satisfies Slater's condition and hence has zero duality gap. 
Furthermore, toward maximizing the Lagrange function, we 
can generalize (fT5T l into the following form: 



max 

>0,PL. 



W k 4>r, 



ps pr 
r>/ mnk mnk 
, "= 



)-Xp 



mnk 



l 'ntnk 



"mnk 



Since the partial derivatives of the above maximization ob- 
jective contains P mnk and PL n j, only in the form of ~ mnk 



and 



we always have Pf r * nk and P r s n * nk as the product of 



4>mnk and a non-negative factor. This leads to a maximization 
problem of the form in ([T8l l. which has been shown to admit 
a binary optimal solution in Section IIV-CI 

Besides DF, the time-sharing variants of any relaying strate- 
gies with long-term or short-term average power constraints, as 
well as all capacity achieving strategies, have concave achiev- 
able rates 0341 , Our algorithm is applicable to these current 

4 Consider the following idealized example for illustration. If = \ for 
all I, the number of iterations would need to be L = 0(e^ mal! ) to satisfy 
the convergence requirement X)z=i = 6(Amai). 



and future relaying strategies to find the optimal solution. 
However, the closed-form solutions for (P^*,,,P"i.) may 
be difficult to find in some cases, requiring more involved 
numerical computation. 

For relaying strategies that do not have concave achievable 
rates, such as AF, near-optimal solutions can be obtained by 
using the proposed approach in the following senses: 

• A concave bound of the achievable rate may be used to 
approximate R(m,n,k). For example, with AF, we have 
R(m, n, k) = i log(l + 1+ ;"; b ;|g;gf3?£ ifc + P mnk cmk)- 
A concave upperbound is obtained by removing "1" from 
the denominator. By substituting such a concave bound 
for R(m, n, k) in the original optimization problem, we 
obtain a solution that optimizes in terms of the bound. In 
the case of AF, such solution is near-optimal for weighted 
sum-rate, since the "1" is negligible for paths with high 
effective SNR, while paths with low effective SNR do not 
contribute substantially to the performance objective. 

• It has been shown in l25ll that, regardless of the convex- 
ity of the objective function in a multi-channel resource 
assignment problem, if the objective at optimum is a 
concave function of the maximum allowed powers, the 
duality gap of the Lagrange dual induced by power con- 
straints is zero. This is due to time-sharing over resource 
assignment strategies. Furthermore, there is a frequency- 
domain approximation of time-sharing, so that the duality 
gap is asymptotically zero when the number of channels 
goes to infinity. Hence, for systems with a large number 
of channels, near-optimal results can be achieved by the 
proposed approach. 

Finally, if we consider R(m, n, k) as a general concave 
utility function of the rate on path "P(m, n, k), then for concave 
and increasing rates, the utility function is also concave 
in the optimization variables, so that a similar optimization 
approach is applicable. An example is the weighted a-fair 
utility function J35), which represents general fairness targets 
such as proportional fairness and max-min fairness. Note that 
this provides only fairness among the paths, instead of among 
users who could be assigned multiple paths. 



VI. Numerical Results 

In this section, we compare the performance of jointly 
optimal channel pairing, channel-user assignment, and power 
allocation with that of suboptimal schemes. We further study 
the different factors that affect the performance gap under 
these schemes, in order to shed light on the tradeoff between 
performance optimality and implementation complexity. The 
suboptimal schemes considered are 

• No Pairing: Channel-user assignment and power allo- 
cation are jointly optimized, but no channel pairing is 
performed, i.e. the same incoming and outgoing channels 
are assumed. The solution is found by always assigning 
an identity matrix to X instead of solving (l24l i. 

• No PA: Allocate power uniformly across all chan- 
nels, subject to power constraints. Channel pairing and 
channel-user assignment are jointly optimized by solving 
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user 1 




user 2 



user 3 



user 4 

Fig. 2. Simulation configuration with K = 4 users 

(O with given power, which is a three-dimensional as- 
signment problem over binary *&. The solution is found 
by following the procedure in Section IIV-CI 
• Separate Optimization: The three-stage solution proposed 
in lfT8l . with channel-user assignment based on maximum 
channel gain over the second-hop channels, channel pair- 
ing based on sorted channel gains, and water-filling power 
allocation. 

« Max Channel Gain: Channel-user assignment by maxi- 
mum channel gain over the second hop, with uniform 
power allocation and no channel pairing. 

We use OFDMA as an example for a multi-channel system. 
The relaying network setup is shown in Fig. [2] where the 
distance between the source and the relay is denoted by d sr , 
and K = 4 users are located on a half-circle arc around the 
relay with radius d r d- A 4-tap frequency-selective propagation 
channel is assumed for each hop, and the number of channels 
is set to N = 16. We define a nominal SNR, denoted by 
SNR nom , as the average received SNR over each subcarrier 
under uniform power allocation. Specifically, with total power 

constraint P t , we have SNR nom = ^o-^jv — ' where n = 3 
denotes the pathloss exponent, a 2 denotes the noise power 
per channel, and d s d denotes the average distance between the 
source and users. A total power constraint and equal individual 
power constraints on both the source and the relay are assumed 
with P s = P r = |Pt, unless it is stated otherwise. 

A. Performance versus Nominal SNR 

We compare the performance of various channel assignment 
and power allocation schemes at different SNR nom levels 
for K = 4. We fix the ratio d sr /d r d to be 1/3. Fig. [3] 
depicts the normalized weighted sum-rate (normalized over 
N) vs. SNR nom for DF relaying with equal weight, i.e., w = 
[wk]ixK = [.25, .25, .25, .25]. The jointly optimal scheme 
outperforms the other suboptimal schemes, and provides as 
much as 20% gain over the Separate Optimization scheme. 
The gain is increased when an unequal weight vector w is 
required to satisfy different user QoS demands or fairness. 
Fig.|4]shows the normalized weighted sum-rate vs. SNR nom for 
w = [.15, .15, .35, .35], where a substantial gain is observed 
by employing the jointly optimal solution. 

B. Performance versus Number of Users 

In this experiment, we show how the number of users affects 
the performance of various resource assginment schemes. We 
increase the number of users in Fig [2] and uniformly place 




Nominal SNR, dB 

Fig. 3. Normalized weighted sum-rate vs. nominal SNR with w = 
[.25, .25, .25, .25], N = 16, K = 4, and DF relaying. 




0.1 1 1 1 1 1 1 1 1 1 

012 3 456789 

Nominal SNR, dB 

Fig. 4. Normalized weighted sum-rate vs. nominal SNR with w = 
[.15, .15, .35, .35], N = 16, K = 4, and DF relaying. 

them around the half-circle arc. In order to properly compare 
weighted sum-rate under different number of users, we do not 
normalize , Wk = 1. Instead, we ixwi = l for all k. The 
nominal SNR is SNR nom = 4dB, and the ratio d sr /d r d = 1/3. 
Fig |5] shows the normalized weighted-sum rate vs. the number 
of users for DF relaying under total power constraint P t . As we 
see, the sum-rate is improved due to the multi-user diversity 
gain with an increased number of users. In addition, consistent 
performance gain under joint optimization can be seen over 
different user population sizes. 

C. Impact of Relay Position 

Through this experiment, we study how the relay position 
affects the performance under various resource assignment 
schemes. The K = 4 users are located close to each other 
as a cluster, and they have approximately the same distance to 
the relay and the source. We change the relay position along 
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- Joint Opt. 

- No Pairing 

- No PA 
Separate Opt. 

- Max Channel Gain 



15 20 
No. of Users 




Fig. 5. 

m = l 



Normalized weighted sum-rate vs. number of users with equal weight R 6 Normalized weighted sum . rate vs . relay location; K = 4 , 
for 1 < i < A, N = 16, and DF relaying. ^ ^ 3^ 3^ N = ^ md Dp relaying 



the path between the source and the user cluster. Figs [6] and 
[7] demonstrate the normalized weighted-sum rate vs. the ratio 
d sr /d sd . We set w = [.15, .15, .35, .35], and SNR nom = 3dB. 
Fig. [6] shows the DF relaying case under both total and 
individual power constraints, and Fig.|7]shows the AF relaying 
case under a total power constraint. 

We see from Fig. [6] that better performance is observed 
when the relay is closer to the source than to the users in 
DF relaying, as correctly decoding data at relay is important 
in successful DF relaying. In addition, comparing the joint 
optimal scheme with No Pairing scheme, we see that the 
gain of channel pairing is evident when the relay is closer 
to the source, but diminishes when the relay moves closer 
to the users. In the latter case, as the first-hop becomes the 
bottleneck, channel pairing at the second-hop provides no 
benefit. This is not the case for AF relaying. As shown in 
Fig. [7] channel pairing gain is observed throughout different 
relay positions. Furthermore, the performance of the jointly 
optimal solution only has mild variation throughout different 
relay positions, unlike the No PA scheme. This suggests that 
the benefit of optimal power allocation for AF relaying is more 
significant when the relay is closer to either the source or the 
users. 

VII. Conclusion 

We have studied the problem of jointly optimizing chan- 
nel pairing, channel-user assignment, and power allocation 
in a general single-relay multi-channel multi-user system. 
Although such joint optimization naturally leads to a mixed- 
integer programming formulation, we show that there is an 
efficient algorithm to find an optimal solution to our problem. 
The proposed approach transforms the original problem into 
a specially structured three-dimensional assignment problem, 
which not only preserves the binary constraints and strong La- 
grange duality, but in some cases can also lead to polynomial- 
time computation complexity through careful choices of the 
optimization trajectory. The proposed framework is applicable 




d Id 

sr Si 



Fig. 7. Normalized weighted sum-rate vs. relay location; K = 4, w = 
[.15, .15, .35, .35], N = 16, and AF relaying. 

to a wide variety of scenarios. The potentially significant im- 
provement of system performance over suboptimal alternatives 
demonstrates the benefit of judicial design in such systems. 

Appendix A 
Derivation of Equation ( fT~6b 

For notational simplicity, we drop all subscripts m, n, and 
k from ( fT~5T >. We have the following maximization problem, 
which can be solved in the two cases below. 

w~ t / aP s \ 
max —0 mm | log H — -~— j , 

/ cP s bP r \ -1 
log 1 + + j - (A s + X t )P s - (A r + \ t )P r 

\ 9 <p J ' 

(29) 

s.t. P s ,P r >0. 
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A. Case One : a < c 

In this case, the first term inside the min function in d29t is 
always smaller than the second term. Hence, (|29l i is reduced 
to 

w ~ ( aP s \ 
max -01og ^1 + — J - (A. + X t )P s - (A r + X t )P r 

s.t. P s ,P r >Q. (30) 
Then, the optimal solutions from water-filling are obtained as 



ps* _ 



1 



2(A S + A t )ln2 a 



pr* = q 



(31) 



B. Case Two : a > c 

For this more complicated case, we propose the following 
solution. We inspect the two possible outcomes in comparing 
the first and second terms in the min function in d29l ) at 
optimality. Two separate maximization of (|29b are performed 
under the constraint of either outcome. Then, the optimal 
(P s ,P r ) is given by the better of these two solutions. 

1) Assumption 1: aP s * < bP r * + cP s * : Under this 
assumption, we have b > and the following optimization 
problem: 

aP s 

(A s + X t )P s - (Ar + X t )P r 



max — 01og ( 1 + 

s.t. (i)P s ,P r >0 

(ii) aP s < bP r + cP s 



(32) 



It has two possible solutions from the KKT conditions. One 
is obtained when the Lagrange multiplier corresponding to 
constraint (ii) is zero and the constraint is strictly satisfied. 
This implies that 

However, this solution contradicts with the assumption that 
(ii) is strictly satisfied. The other, correct solution occurs at 
the border P r = ^^P s . By inserting this into the objective 
function, we have 



2b(X s + A t ) In 2 + 2(a - c)(X r + A t ) In 2 a 



(34) 



In the former case, b > since a ^ c. Furthermore, since 
the first and second terms in the min function in (|29| | are the 
same, we obtain the same solution as in (f34t . 

In the latter case, we define two new variables V s — (A s + 
X t )P s and V r = (X r + X t )P r . Substituting them into the 
objective of (I35K we have 

id- / cV s bV r 

max — log 1 



y, yr 2 



0(A S + A t ) 0(A r + A t ) 



- V s - V r . 



The solution depends on the relation between 
b . 

. If 



(36) 
and 



> 



, then we have V r * = 0, since 
otherwise a better solution to (1361 1 would be (V s = 
V s* + yr*^ V r = Q y Substituting V r * = into (O, 

_ !L _ _ A a +A t 1+ 7 
2 In 2 c J r- 

If A7TI7 = a^aI' (® is a function of (V s + V r ) only, 
and V r * = is a maximizer. Hence, again we have 

t/s* r w _ x s +\ t i ^ ~ 

v ~ L 2 In 2 
If 



we have V s * = 



< 



A -, similarly we have V s * = 0. However, 
this together with our assumption that constraint (ii) of 
( 1351 1 is strictly satisfied, i.e., bP r * < (a — c)P s * , implies 
that V r * < 0, which is not a feasible solution. Therefore, 
in this case at optimality the condition bP r * = (a—c)P s * 
prevails. 

Appendix B 
Proof of LemmaQ] 

Given any * = [4> mn k]NxNxK with < 4> mnk < 



1 and satisfying ([TT1 >. let 



E 



From 



m— 1 Z_^fc= 



we have Y^, 
can be constructed as 



= 1, we have J2 m =i x mn = 1. Similarly, 
= 1. Hence < x mn < 1. Then, y% n 



ran 
Vk 



1/K, 



> o 
= o 



(37) 



P' 



Hence, EfeLi^" = 1 and < yf n < 1. Note that 1/K 
above is arbitrarily chosen, and the mapping from <& to (X, 
y mn ) is one-to-many. 

Given X and y mn with < x mn < 1 and < y^ n < 
1, satisfying E ra =i x mn = l,Vm, E^i^mn = ljVn, and 
EfcLiy™" = l,Vm,n, clearly < ^ mrifc = x mn y£ n < 1, 
and it is easy to verify that ( fTTT i is satisfied. This establishes 
the equivalence of 4> and the proposed decomposition. 



2) Assumption 2: aP s * > bP r * + cP s * : Under this 
assumption, we have the following optimization problem: 

w ~ ( cP s bP r 
max — log 1 H = — I — 

p^p- 2 r ^ 

(A, + X t )P s - (X r + X t )P r 
s.t. (i)P s ,P r >0 

(ii) aP s > bP r + cP s (35) 

From the KKT conditions, at optimality, either bP r * = (a — 
c)P s *, or the Lagrange multiplier corresponding to constraint 
(ii) is zero and the constraint is strictly satisfied. 



Appendix C 
Proof of Lemma|2] 

We note that there exists at least one index vector 
(m', n', k') such that 0* n < n /fc'(A*) = 1. Furthermore, this cho- 
sen path must have non-degenerate user weight and channel 
gains so that the weighted rate function Wk>R(m' , n', k') is not 
uniformly zero, i.e., Ww > 0, a m ' > 0, and b n >k> +Cm<k' > 0. 

Suppose there exists 4>m'rik'(^*^ = 1 such that 



P% n >k>W is either [< 



i - 



or 



— — . In the former < c, n 'k>, and the latter, 

a-' and > x"+x* • Furthermore, since b n 'k> 
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and Cm'k' cannot both be zero in the latter case, c m i k > > 0. 
Then we have 

min{P s ,P t } > £ f£* nfc (A*) > P^ W (A*) 



> min 



iin{[ 



w k > 



1 

a{\* +X* t ) _ 
ww 1 -1+1 

a(X* + X*) c m , fc ' J J 
w k > 1 



otherwise the achieved sum-rate is uniformly zero. Then 

for any (rn',n',k') such that §* m i n i k , (A*) = 1, 

min{P s ,Pj > 2 P°; nk (\*) > i£W**) 
w k >b n > k > 1 



> 



a(6 n 'fc'(A* + A*) + a m < (A* + A t *)) a, ; 



(43) 



> 



a(A*+A^) mm{a m ,,c m >k>} 



(38) 



Hence, A? + A* > . "V min{a ,,c ,»,} A * ^ 

' s t — ct(mm{a m / ,c m i k i } min{P s ,P t } + l) 1 

Otherwise, for all ^ nfc (A*) = 1, we have a m > c mk and 



DS* (\*\ - WkOnk _ _ 

proceed with the following cases: 



i 1 + 



We 



If there exists 



< 



rnin{P s ,P t } > 



> 



6* n /„/ fc /(A*) = 1 such that c m < k - > and 
then we have 



-f?rmfc(A ) > P^, n , k ,(X ) 
W k /b n ' k ' 



> 



a(o«'fc'(A* + Aj) + (aw - c m > k ')(X* + A|)) 

Wk'b n 'k' 

b n > k '( a m> 



which implies that (A* + A,*) + jrf^iX* + K) ^ 
, . r n Wh ' : ; — r^r. Considering the extreme case for the 

a(mm{Ps ,-Pt }+ - — ) 

slope and intercept of this linear inequality, we have 

A* e n 2 . 

Appendix D 
Proof of Lemma[3] 

To find an upper bound for |j A*||2, we consider the follow- 
ing cases for the activation pattern of the individual power 
constraints (0 and the total power constraint (@) at global 
optimum. 

• Neither constraint in (01 is active: In this case, (0 must 
be active, since otherwise there would be more power to 
increase the sum-rate. Thus, we have X* = A* = and 



1 



Pt = 



(W(A: + A*) + °--«yy-y«> {K + A|)) a m . 



m,n,k 



Pmnki^ ) + Pmnki^ ) 



(44) 



(39) 
, and 



which implies that A! + A? > —> h ' f '"' fc ' , , ^ 

hence A* e TL\. 

Else, if there exists (l>m'n'k' (A*) = 1 sucn that 

x?+x* > x"+\" 1 then Cm ' k ' > since bn ' k ' and 
c m 'fc< cannot both be zero. In this case, the third ex- 
pression in (fTol l applies to the path (m! ,n' , k'). Since 

[ a (b n , kl {\t+\;)+\*2> k -c m >k>KK+K)) ~ v^"] 1S the op " 
timal power allocation for this path, we have 

r , ,, P el p rl \*\ ~> 

'-'m'n'k' A 1 ? r m'n'k' ' £ m'n'k' ) ^ I — 

£m'rik' (lj P m 'rik' ! Pm'n'k' i A ) , (40) 



Since C„ fc (A*) < 1, substituting = 1 and 

A* = A* = into ( fl6b . and considering all possible 
scenarios of ( fTol l. we have 



Pt < J]] max{ [ 



Wk_ J- 

aAJ a„ 



1 



} 



E 



w fc _ A 2 
aA? aXt 



L aA( c mfc J ' aAj 
(45) 



Hence, we have A? < so that 



< 



N 2 
aP 



(46) 



where 



m'n'k' 



ps2 a 

m'n'k' 



Wk'b n > 



a(bn' k '(X* s + Af ) + (aw - <W fe ')(A;!: + A?)) 
1 



Bof/i constraints in (0 are active: We have 



niAi.k 



Pr — 22 Pmnki^* 



(47) 



1 



rn,7i,k 



"(As + A?) c m ' fe ' 



(41) 



correspond to the cases a m 'Pm'n'k' 



pri 



m'n'k' u n' k' 1 m 'n'k' 

c m'k'Prn'n'k'' respectively. This implies that (j42j>. 
Hence, 



bn'k' PL 2 



Again, substituting ^^^(A*) = 1 into ( TTBT l. and consid- 
ering all possible scenarios, we conclude that 

1 1 + 



P s < max 



a: + x* t > 



Wk'C m ' k ' 



m,n,k 
Wk 



-a{X* s +\* t ) a Tl 

1 n + 



4a(a m -min{P s ,P t } + l)' 



a(A* + A t *) Cmk 1 ' 

l(a m > c mk )w k b nk 



so A* e ^i. 

Else, the only scenario left is when c mk = for {m, fc : 
4>mnk(^*) = in m is case, there exists b nk > 0, since 



J))} 



< 



a(MA* + A t *) + {a m - c mfc )(A* + A*)) 
A 2 



E 



a(A* + A t *) a(A*+A t *) ' 



(48) 
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^ log(l + a m ,P£ n , k ,) (K + A*)(l + am \ Cm ' k ' )P^, n . k . > log(l + c m , k ,P% n , k ,) (A 5 * + X* t )P^ n , k , 

A O n 'fe' A 

^ log(l + a m .P% n , k .) + (\* s + \* t )P%, n . k , > ^ log(l + c m , k ,P:3 n , k ,) 

log(l + a m , min{P s , P t }) + ^L>^E_ i og (l + c m , k ,P^ n , k ,) 
A a A 

4 + 4<wmin{P s ,P t } > 1 + 

4 + 4<w min{P , P t } > ^"TTm " (42) 



P r < 



l(a m > c mk )w k (a m - c mk ) 



^, k a ( b ^(K + K) + (am - Cmk)(X* + A*)) 

< y Wk = N2 (49) 



Hence, we have A! + Xt < -V, and A* + AJ < so 
that 

V2N 2 

l|A*|| 2 < V j p V p , ■ (50) 
amm{P s , P r | 

• Only one constraint in (0) is active: We have either 
A* = and (08), or A* = and d49l . For both cases, an 
upperbound for ||A*||2 is given by d50t . 
Summarizing the three cases above, we have 

V2N 2 



< 



amin{P s ,P r ,P t } 



(51) 



Appendix E 
Proof of Lemma|4] 

We first note that there must exist one path (m, n, k) with 
non-degenerate user weight and channel gains so that the 
weighted rate function w k R(m,n,k) is not uniformly zero, 
i.e., w k > 0, a m > 0, and b nk + c mk > 0. In Algorithm Q] 
subgradient updating is performed over two cases, either there 
exists some c mk > and X 1 - 1 ^ G IZi, or c mk = for all m 
and k and G K 2 . 

In the former case, let t\ = 
min Wk min{ min a m , min c mk } 

{fc:to fc >0} {m:a m >0} {m, k :c mk >0} ™ 

4a(maxa m min{P s , P t } + 1) 

m 

projected subgradient updating is performed over 
As + Xf- > £i. Since there exists some m and k 
such that Wfe > 0, a m > 0, and c mk > 0, we have ei > 0. 
From (1161 . we see that in all scenarios 



rnnk 



(A«) < 



< < oo 



^r„ fc (A (i) ) < 



p „ (A(I)) < Ma m -c mfc ) 

ae\b nk 



The second inequality above hold since Pmnki^ 1 ) = when 

bnk = 0. 

In the latter case, let e 2 = 
min Wk 

{k:w k >Q\ 

a( min{P s , P t } + : 

mm a„ 

{m:a m >0} 



(52) 



< oo 



(53) 



The projected subgradient 



mm a m 

updating is performed over (A a + A t ) + {m ^\ (A r + 

n,k 

A*) > £2- Since there exists some m and k such that wj. > 0, 
a m > 0, we have e 2 > 0. Furthermore, b n k > since From 
(TToT l. we see that in all scenarios 



^*„ fc (A (i) ) < 



»7v 



a(xi l) + f-A^) «e 2 



< < oo , 



- r mnfc( A I ^ 7 - r mnfc( A I 

Onk 



ajbnk 



(54) 



< OO 

(55) 



The second inequality above hold since P m * n k^^) = when 

fonfc = 0. 

We note that the bounds in (f52T>-d55T> are not functions of 
m and n. Substituting the bounds of either of these two cases 
into d27t , and noting that only N 2 paths are chosen, we have 

||0(A)|| 2 = O(7V 2 ). 
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